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ABSTRACT 

The  newsboy  problem  solutions  under  the  conditions  of 
risk  and  of  uncertainty  about  demand  are  well  known.   The 
former  is  the  case  where  the  distribution  of  demand  is  known 
or  estimated,  and  the  latter  is  the  case  where  it  is  not 
known  but  the  range  of  demand  is  given. 

Minimizing  the  expected  cost  under  risk  and  the  minimax 
approach  under  uncertainty  are  well-known  methods  to  solve 
the  problem.   The  situation  is  considered  where  demand  fre- 
quency data  is  acquired,  one  observation  per  decision  period. 
For  the  expected  value  solution,  this  study  presents  five 
candidate  estimators  based  on  the  order  statistics  and 
evaluates  them  by  comparing  costs  with  those  achieved  using 
the  minimax  rule  under  uncertainty.   This  is  done  by  simu- 
lation; the  results  also  provide  information  about  the  switch- 
ing period  from  the  minimax  rule  to  the  proposed  quantile 
estimators. 
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I.   INTRODUCTION 

The  newsboy  problem  was  described  by  Morse  and  Kimball 
early  in  1950  [1].   This  problem  has  been  presented  in  the 
literature  under  a  variety  of  names,  including  newsvendor 
problem,  Christmas  tree  problem,  and  high-fashion  apparel 
problem.   This  kind  of  problem  may  be  described  in  the 
following  way:   The  newsboy  has  to  decide  how  many  papers  to 
purchase  for  resale.   If  he  buys  too  many  papers  he  will 
incur  loss  due  to  left-over  papers  (little  salvage  cost) . 
If,  on  the  other  hand,  he  buys  too  few,  he  should  also  incur 
an  opportunity  cost  due  to  lost  sales.   To  determine  the 
optimal  order  quantity,  it  is  necessary  to  balance  two 
opposing  costs.   Solutions  to  the  newsboy  problem  under  both 
conditions  of  risk  and  uncertainty  are  well  known.   The 
determination  or  estimation  of  the  demand  distribution  is 
critical  in  dealing  with  this  problem.   Optimal  solutions  to 
the  newsboy  problem  which  minimize  expected  cost  require 
information  about  a  particular  quantile  of  the  demand  dis- 
tributions.  If  the  demand  follows  some  known  distribution 
the  problem  can  be  solved  by  decision  making  under  risk,  as 
is  also  the  case  when  the  key  quantile  is  known.   With  only 
information  concerning  the  range  of  demand,  optimal  solutions 
can  be  found  by  some  procedure  such  as  the  minimax  rule  under 
uncertainty. 

The  purpose  of  this  study  may  be  summarized  as  follows. 
We  are  concerned  with  decision  policies  during  a  period  when 


demand  data  is  being  acquired  at  the  rate  of  one  observation 


per  decision  period,  starting  with  no  data  except  an  estimate 
of  the  maximiom  demand.   To  find  the  optimal  expected  value 
decision  rule,  five  candidate  quantile  estimators  will  be 
proposed  on  the  basis  of  order  statistics.   They  will  be 
evaluated  by  comparing  their  performance  with  each  other  and 
with  the  well-known  minimax  rule  under  uncertainty.   An  esti- 
mator which  has  the  smallest  relative  cost  deviation  to  the 
risk  cost  solution  (ideal  case)  will  be  proposed  as  optimal. 
We  shall  also  examine  the  switching  time  to  change  from  the 
uncertainty  rule  to  applying  the  proposed  estimator. 

In  Chapter  II,  we  shall  review  the  newsboy  problem  solu- 
tions under  conditions  of  both  risk  and  uncertainty,  and  then 
suggest  alternative  approaches  for  demand  estimation,  and 
for  the  switching  period  from  uncertainty  to  risk.   For  esti- 
mation, five  quantile  estimators  will  be  proposed  in  Chapter 
III.   In  Chapter  IV,  candidate  quantile  estimators  will  be 
evaluated  by  simulation.   The  performance  of  the  estimators, 
and  candidate  switching  periods  will  be  discussed  in  Chapter 
V.   Conclusions  and  recommendations,  for  further  work,  will  be 
given  in  the  final  chapter. 


II.   STATEMENT  OF  THE  PROBLEM 

In  this  section  we  shall  briefly  describe  the  well-known 
newsboy  problem,  illustrating  its  optimal  solutions  under 
both  risk  and  uncertainty  with  respect  to  demand.   Then  we 
shall  describe  the  quantile  estimation  problem  we  wish  to 
solve. 

A.   THE  NEWSBOY  PROBLEM 

The  problem  may  be  generalized  in  the  following  way. 
The  decision  maker  selects  a  number  (the  newsboy's  order 
quantity)  but  the  actual  demand  (the  number  of  papers  his 
customers  will  want  to  purchase)  may  have  a  randomly  deter- 
mined value.   If  the  decision  maker's  number  is  greater  than 
this  he  pays  a  cost  proportional  to  the  difference,  or  he 
pays  a  cost  proportional  to  the  amount  by  which  his  number  is 
smaller  than  this.   His  cost  is  zero  only  if  his  number  is 
the  same  as  the  actual  demand. 

Let  demand  D  be  a  continuous  random  variable  and  let 
quantity  S  be  the  number  the  decision  maker  selects  to  have 
on  hand.   The  newsboy  cost  equation  may  be  formulated  as  a 
two-piece  continuous  linear  function  in  the  following  manner. 


C(S)   = 


I 


C  (S  -  D)     if     0  <_  D  <^  S 


(1) 
^  Cq(D  -  S)     if     D  >  S, 


where  C   is  unit  cost  of  surplus  when  S  exceeds  D,  and  C 

s  '^  o 

is  unit  cost  of  outage  or  shortage  when  D  exceeds  S.   The 
cost  function  is  shown  in  Figure  1. 


Demand^  D 


Figure  1.   The  Newsboy  Cost  Function 

B.   THE  NEWSBOY  PROBLEM  SOLUTIONS  UNDER  RISK 

For  our  case  where  demand  is  a  continuous  random  varia- 
ble D  with  probability  density  function  f (D) ,  and  S  is  a 
decision  variable  whose  value  the  decision  maker  is  trying 
to  find  over  D,  the  expected  cost  may  be  represented  as: 


E[C(S)]   =   C    I  (S  -  D)f  (D)dD  +  C    f  (D-S)f(D)dD. 
s   ^  o   ^ 

0  s 

(2) 


We  wish  to  find  quantity  S,  which  minimizes  the  expected 
cost.   The  optimum  (minimum  in  this  case)  value  of  S  may 
be  found  by  differentiating  the  expected  cost  function  (2) 
with  respect  to  S,  setting  the  derivative  to  zero,  and  solving 
for  the  optimum  value  S*: 


"^^^^^S^^^^   =   C^   /  f  (D)  dD  -  C^   I  f  (D)  dD   =   0.  (3) 

0  s 


From  this  we  obtain  the  optimal  value  S*  using  the  cumula- 
tive distribution  function 


X 
F(X)   =    /  f(D)  dD, 
0 


and  the  result  is  the  well-known  result. 


s   o 


Differentiating  Equation  (3)  with  respect  to  S  to  get  the 
second  derivative  at  S* , 


d^E[C(S)] 

d-s2 


=   (C^  +  C  )  f(S*)  . 
s=s*       ^   ° 


Hence,  S*  is  the  minimum  expected  cost  solution  since 

(C  +C  )f(S*)  is  positive. 
s    o 

To  summarize,  the  optimal  solution  S*  to  minimize  the 
expected  cost  is  the  [C  /(C  +  C  ) ] th  quantile  of  the  demand 
distribution,  as  shown  in  Figure  2,  and  the  crucial  infor- 
mation we  need  about  the  distribution  is  that  quantile. 


C.   THE  NEWSBOY  PROBLEM  SOLUTIONS  UNDER  UNCERTAINTY 

The  decision  maker  may  be  confronted  with  a  newsboy- 
problem  without  knowing  or  being  able  to  estimate  the 
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q  = 


Demand ,  D 


Figure  2.   Minimum  Expected  Cost  Solution 


probability  distribution  of  demand.   The  problem  of  deciding 
on  S  in  this  situation  is  called  a  decision  under  uncertainty 
The  decision  maker  may  be  able  to  estimate  an  upper  bound 
for  the  demand  value.   In  this  case  we  say  that  all  we  know 


\ 


about  demand  is  that  demand  D  is  such  that  0  <  D  <  D 

—   —  max 

S^  Three  well  known  approaches  to  decision  making  under 
uncertainty  are  Laplace  solution  under  assumption  of  a  uni- 
form distribution  of  demand,  minimax  cost  solutions  and 
\  "'minimax  regret  solutions  [2]  .   For  the  newsboy  problem, 

these  all  lead  to  the  same  rule  to  find  the  optimal  solution. 
We  shall  discuss  only  the  minimax  cost  approach.   This 
approach  is  to  choose  S  so  that  the  worst  possible  cost  will 
be  minimized,  that  is,  to  minimize  maximum  cost.   The  maxi- 
mum cost  occurs  at  D  =  0  or  at  D 

cost  coefficients.   This  can  be  seen  from  Figure  1.   If 
maximum  cost  occurs  at  D  =  0 ,  we  would  want  to  reduce  S 


D    ,  depending  on  the 
max    ^     ^ 
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while,  if  maximum  cost  occurs  at  D  =  D    ,  we  would  want 

max 

to  increase  S  to  reduce  cost.   We  can  minimize  the  maximum 

cost  when  S  is  chosen  so  that  the  cost  at  D  =  0  is  equal  to 

the  cost  at  D  =  D    .   Equating  these,  we  obtain: 

max 


C   S*   =   C  (D    -  S*) 
s         o  max 


or 


S*  =   'ra-'  °max  •    ^"  <  =  ' 


So  far,  we  have  introduced  the  minimum  expected  cost 

solution  under  both  risk  and  uncertainty  conditions.   The 

minimax  approach  is  tied  to  our  initial  estimate  of  D    , 
^^  max 

and  thus  the  optimal  value  S*  from  the  minimax  rule  will       ^^4A 


; 


be  in  error  to  the  extent  that  our  estimation  of  D    is 

max 

in  error.   This  error  will  continue  over  many  decision 

periods  unless  our  estimate  of  D    is  revised. 
^  max 

D.   QUANTILE  ESTIMATION  AND  A  DECISION  RULE  SWITCHING 
POINT  FOR  THE  NEWSBOY  PROBLEM 

We  are  interested  in  finding  a  quantile  estimator  which 
gives  better  results  than  the  minimax  approach  under  condi- 
tions of  limited  demand  data,  and  also  in  finding  the 
switching  period  which  is  the  time  at  which  we  switch  from 
the  minimax  rule  to  the  expected  cost  rule  using  the  quantile 
estimator.   The  purposes  of  this  study  are  to  find  the  "best" 
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quantile  estimator  among  five  candidate  estimators,  and  to 
obtain  information  about  the  switching  period. 

We  recall  that  we  wish  to  find  the  quantile  of  the 
cumulative  distribution  function,  F(D).   The  optimal  S*  is 


found  by  an  inverse  mapping  of  the  [C  /(C  +C  )  ]  th  quantile 

of  F(D)  into  demand  D  under  risk  condition.   Let  q  =  C  / 

^    o 

(C  +  C  ) ,  then  we  have: 
s   o 


S*   =   F"^(q)  . 


Although  we  don't  know  the  demand  distribution,  each  suc- 
cessive decision  period  provides  another  observation  from  the 
demand  distribution,  and  we  will  assume  that  an  observation 
of  demand  is  indeed  obtained  each  period.   As  we  acquire  more 
data  we  might  be  able  to  estimate  S*  =  F   (q)  as  a  function 
of  the  given  data.   Two  kinds  of  approachs  are: 

1.  a  parameteric  approach  to  hypothesizing  the  form  of 
demand  distribution, 

2.  a  non-parametric  approach  to  find  the  quantile. 

For  the  first  case,  the  normal,  and  exponential  distribu- 
tions have  been  found  to  be  of  considerable  value  in  des- 
cribing demand  distributions.   The  normal  distribution  has 
been  found  to  describe  many  demand  functions  at  the  factory 
level;  the  exponential,  at  the  wholesale  and  retail  levels 
[3] .   Of  course,  these  distributions  should  not  be  auto- 
matically applied  to  any  demand  distribution.   Statistical 
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tests  should  establish  the  basis  for  any  standard  assumption 
concerning  a  demand  function. 

We  are  directing  our  attention  to  the  case  where  the 
expected  value  solution  is  not  applicable  because  of  lack  of 
available  information.   Further,  if  we  had  enough  information 
for  a  fit  of  a  distribution,  we  should  have  a  good  estimate 
of  the  quantile  we  seek.   Hence,  the  second  case,  the  non- 
parametric  approach,  will  be  considered  to  find  the  optimal 
S*  and  switching  period. 
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III.   DEVELOPMENT  OF  THE  ESTIMATORS  ON 
THE  BASIS  OF  ORDER  STATISTICS 

As  mentioned  in  the  previous  chapter,  we  are  subject 
to  find  the  optimum  (minimum  cost  solution)  of  a  newsboy 
problem  without  knowing  the  q-th  quantile  of  the  demand  dis- 
tribution.  Consequently  a  non-parametric  or  distribution- 
free  approach  is  suggested  to  find  the  solution.   In  this 
chapter,  we  will  introduce  the  notion  of  order  statistics  as 
a  basis  for  estimating  this  quantile,  and  suggest  five  can- 
didate quantile  estimators  which  may  be  used  in  newsboy 
solutions . 

A.   ESTIMATORS  USING  ORDER  STATISTICS 

The  non-parametric  approach  we  shall  follow  is  on  the 
basis  of  order  statistics  and  although  it  is  applicalbe  to 
both  continuous  and  discrete  random  variables,  we  shall  direct 
our  attention  to  the  continuous  case.   We  are  interested  in 
the  case  where  the  population  is  not  known,  but  we  have  a 
random  sample  of  size  n  from  its  unknown  cumulative  distribu- 
tion function  F(X).   Let  X, ,X_,...,X  be  a  random  sample. 

1   2       n  ^ 

Then  X(l)  <_  X(2)  1  ...  ^  X(n)  ,  where  the  X(i)  are  arranged 
in  order  of  increasing  magnitudes,  and  are  defined  to  be  the 
order  statistics.   Let  X  be  the  q-th  quantile  of  the  popu- 
lation.  The  estimate  of  X   can  be  expressed  as  a  function  of 


I '  9 / • • • / X  ,  I.e., 
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Estimate  of  X   =   f (X^ ,X^ , . . . ,X  )  . 
q        12      n 


We  shall  consider  now  a  function  F  (x)  which  is  construc- 

n 

ted  from  the  sample  values  of  X, ,X»,...,X  .   The  function 
'^  1   2      n 

F  (x)  is  called  empirical  cumulative  distribution  function, 
n 

Let  X, ,x^/...,x   denote  the  observed  values  of  X, ,X^,.../X  . 
12n  12n 

Then  for  each  number  xi-^^    <    x   <    °°)  ,    the  value  of  F  (x)  is 
defined  to  be  the  proportion  of  observed  values  in  the  sample 
which  are  less  than  or  equal  to  x.   It  follows  from  the 
Glivenko-Cantelli  Theorem  that,  as  n  approaches  infinity  the 
empirical  cumulative  distribution  function  F  (x)  converges 
uniformly  to  F(x).   Since  we  are  concerned  with  an  unknown 
cumulative  distribution  function  F(x)  which  is  continuous, 
some  type  of  smoothed  curve  of  F  (x)  might  yield  a  reasonable 
estimator  of  F(x).   Let  q  be  the  [C  /(C  + C  ) ] th  quantile, 
where  we  wish  to  find  the  estimate  of  X   so  that 

q 

F (estimate  of  X  )  =  q.   This  is  represented  in  Figure  3. 
Now,  we  are  interested  in  the  estimator  which  gives  us  mini- 
mum cost,  but  we  have  no  information  about  it  yet.   Let  us 
turn  our  attention  to  the  expected  cost  as  the  function  of  the 
estimate  of  X  .   Let  perturbation  be  the  difference  between 
the  estimate  of  X  and  true  value  of  X  .   Then  the  expected 

q  q 

cost  of  perturbation  is  proportional  to  the  variance  of  the 
given  unbiased  estimator  [4].   The  properties  of  unbiasedness 


16 


cr  = 


Smoothed  F  (x) 

-   F  (x) 
n 


estimate  of  X 


^(r-l)   ^(r)   ^(r+1) 


Demand,  X 


Figure  3.  Estimate  of  X  on  the  Smoothed 
Curve  of  Empi?ical  CDF,  F  (x) 


and  minimum  variance  should  be  considered  for  selecting 
estimators. 

B.   RATIONALE  FOR  THE  ESTIMATORS 

In  this  section,  we  shall  suggest  five  different  quantile 
estim.ators  on  the  basis  of  order  statistics.   In  a  statisti- 
cal evaluation  of  these  estimators,  it  is  desirable  to  use 
the  uniform  (0,1)  distribution,  since  the  expectation  and 
variance  of  its  order  statistics  are  simple  and  exact. 
Investigating  these  properties  of  order  statistics  using  other 
distributions  is  not  easy..  If  we  have  a  distribution  similar 
to  the  uniform,  or  if  we  have  a  large  number  of  observations, 
some  approximations  may  be  possible.   Three  points  which  are 
near  each  other  can  be  assumed  to  be  linear  and  some  estimator 
using  order  statistics  would  be  approximately  unbiased. 
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We  shall  examine  the  properties  of  the  order  statistics 

from  the  uniform  (0,1)  distribution,  and  we  shall  consider 

three  neighboring  order  statistics.   Let  X,  ^.,    X,  .  ,  and 

^  (r-1)    (r) 

^(r+1)  ^^  ^^®  (r-l)th,  r-th,  and  (r+l)th  order  statistics 
respectively.   Define  the  value  r  as: 

r   =   [nq  +  0.5]  ,  (6) 


where  [X]  denotes  the  largest  integer  of  X.   Let  us  consider 
an  estimate  of  X   as  the  linear  combination  of  those; 

q 


(estimate  of  X^)   =   aX^^_^^  +  bX^^^  +  cX^^_^^  , 


where  a,  b,  and  c  are  nonnegative  and  sum  to  unity.   There 
can  be  many  possible  ways  to  take  estimators,  but  we  shall 
examine  typical  ones  considering  the  unbiasedness  and  the 
magnitude  of  variance. 

1.   The  First  Candidate  Estimator;   X 

q 

First  of  all,  the  estimator  should  be  unbiased, 
i.e.,  E[X  ]  =  X  .   If  we  take  the  r-th  order  statistic, 

q    q 

X    as  the  first  candidate  estimator,  i.e.. 


Xg  =   X(^,  ,  (7) 


this  is  an  unbiased  estimator  since  E[X.   ]  =  r/(n+l)  for 
the  r-th  order  statistic.   The  coefficients  are  a  =  c  =  0  and 
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b  =  1  in  this  case.   The  first  candidate  estimator  is  repre- 
sented in  Figure  4. 


F(x) 


1.0 


q  = 


C  +C 
s   o 


Figure  4 


' 

) 

_,--—- 

^^ 

^    q 

y     i 

y 

y 

*^^                              \                   < 

^(r-l)  ^(r) 


X 


Demand^  X 


(r+l) 
The  First  Candidate  Estimator,  X 


2.   The  Second  and  Third  Candidate  Estimators:   X  ,  X 
— — — — — q    q 

We  shall  propose  (r-l)th  and  (r+l)th  order  statistics 
as  the  second  and  third  candidate  estimators: 


^q  -   ^(r-i)  ' 


(8) 


and 


^q  "   ^(r+D- 


(9) 


Even  though  these  are  not  unbiased  estimators,  we  are  not 
certain  which  one  gives  us  better  performance  since  we  are 
considering  the  properties  on  the  uniform  distribution  and 
there  is  some  round  up  error  for  taking  the  value  r. 
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3.   The  Fourth  Candidate  Estimator:   X 

q 

We  shall  propose  an  equally  weighted  linear  combina- 
tion of  the  (r-l)th  and  (r+l)th  order  statistics  as  the 
fourth  candidate  estimator: 


^q  =  l'=^(r-l)  *  ^r+l)i-  ("> 


This  estimator  is  an  unbiased  estimator,  since 


E[X„1   =   i[§^+  §ii)   =    "^ 


q^      2  "-n+l    n+l-"      n+1  ' 
and  has  the  minimum  variance  using 

(^^   ^^)       (n+l)2(n+2) 

of  the  uniform  distribution  [5] .   The  variance  of  this  esti- 
mator is  given  as  follows: 


^(r-1)  +  ^(r+1) ,   _   (n-r) (2r-l)+(r-l) 


Var[X  ]   =   Var[-  ^"  ^'    ' '-±-±^]      = 

^  ^  2(n+l)^(n+2) 


This  estimator  is  represented  in  Figure  5. 

4.   The  Fifth  Candidate  Estimator:   X 

q 

We  shall  propose  an  equally  weighted  linear  combina- 
tion of  the  (r-l)th,  r-th,  and  (r+l)th  order  statistics  as 
the  fifth  candidate  estimator,  or 
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q  = 


^(r-l)  ^(r) 


^(r+1) 


Demand,  X 


Figure  5.   The  Fourth  Candidate  Estimator,  X 


X 


^[X  .   ^  ^  +  X,  ^  +  X  .   ^  J  . 
3   (r-1)     (r)     (r+1) 


(11) 


This  is  also  an  unbiased  estimator  and  the  variance  of  this 
estimator  is  given  by: 


Var  [X  ] 

q 


(n-r) (9r-4)  +  (5r-4) 
9(n+l) ^(n+2) 


The  difference  in  variance  between  the  fourth  estimator  and 

the  fifth  estimator  is  equal  to  t-cTt — ,  t  \  / — r^^^-      When  n  =  1  and 

^       18 (n+1) (n+2) 

-4 
n  =  10  this  difference  is  approximately  0.09  and  4x10   , 

respectively.   The  difference  converges  rapidly  to  zero  as  n 
increases.   Since  the  quantity  difference  may  be  assumed  to 
be  negligible  and  the  variance  is  only  concerned  with  uni- 
form distribution,  it  is  hard  to  say  which  estimator  gives 
us  better  performance. 
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In  summary,  although  there  may  be  many  possible 
candidate  estimators  depending  upon  different  forms  of 
linear  combination,  it  may  be  reasonable  to  examine  five 
previously  proposed  estimators.   They  can  be  shown  with 
properties  as: 


1)  The  First  Estimator:   X   =  X,^.       (unbiased) 

q       \^) 

2)  The  Second  Estimator:   X   =   x;   ,,   (biased) 

q       (r-1) 

3)  The  Third  Estimator:   X   =   X,    ',   (biased) 

q       (r+1) 

4)  The  Fourth  Estimator:   X   =   i[X   _  .  "'■^(''+1)5 

(unbiased  min.  variance) 

5)  The  Fifth  Estimator:   X   =   ^[X,  -,  .  +  X ,    ,+X,   ,J 

q      3'  (r-1)    (r)    (r+1) 

(unbiased) . 


Having  proposed  these  quantile  estimators,  it  is  of 
interest  to  see  how  well  they  perform  in  the  context  of 
newsboy  decision  making,  under  cost  criteria.   This  is  the 
subject  of  the  next  chapter. 
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IV.   EVALUATION  OF  THE  CANDIDATE  ESTIMATORS 

In  this  chapter,  the  five  candidate  estimators  will  be 
tested  and  investigated  for  their  performance  in  the  newsboy 
problem  for  several  demand  distributions.   As  performance 
measures,  two  kinds  of  costs,  namely  an  average  cumulative 
cost  and  a  mixed  cost  which  we  shall  introduce  later  on,  will 
be  obtained  during  the  simulation  and  then  the  performance 
of  the  candidate  estimator  will  be  evaluated  using  the  rela- 
tive cost  deviation  to  the  ideal  cost. 

A.   DESIGN  FOR  SIMULATION 

A  newsboy  simulation  will  consist  of  50  replications, 
each  of  which  simulates  50  decision  periods  for  a  given 
demand  distribution.   The  quantile  value  (q  =  C  /  (C  +C  )) 
for  each  distribution  is  varied  from  0.1  to  0.9  in  increments 
of  0.2. 

To  find  the  value  r  and  consequently  to  find  the  estimate 
of  X  ,  some  modification  of  the  method  of  determining  r  is 
necessary  during  the  first  and  the  last  few  periods,  since 
there  is  an  insufficient  amount  of  data  to  apply  the  candi- 
date estimators.   For  example,  if  the  value  r  for  finding 
the  estimate  is  equal  to  zero,  then  X,,.  will  be  chosen  as 
the  estimate  of  X  ,  and  if  r  =  n,  then  the  estimate  of  X 

q  q 

will  be  X,  >  and  so  forth, 
(n) 

For  each  period,  the  average  cumulative  cost  and  the 
mixed  cost  are  computed  for  seven  different  cases  which  are 
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represented  as  follows: 

1)  Risk  case  (Ideal  solution), 

2)  Uncertainty  case  using  the  minimax  rule, 

3)  A  case  using  the  estimator  X  , 

q 

4)  A  case  using  the  estimator  X  , 

5)  A  case  using  the  estimator  X  , 

q 

6)  A  case  using  the  estimator  X  ,  and 

q 

7)  A  case  using  the  estimator  X  . 

q 

The  cost  for  each  period  is  averaged  over  50  replications,  and 
then  the  average  cxomulative  cost  at  period  k  is  computed  by 
cumulating  the  average  cost  from  the  first  period  to  k-th 
period.   The  mixed  cost  at  k-th  period  is  given  by  summing 
the  average  cost  using  the  minimax  rule  up  to  (k-l)th  period 
and  the  average  cost  using  an  estimator  from  the  k-th  period 
to  the  last  period. 

B.   DISTRIBUTION  AND  PARAMETERS 

As  mentioned  previously,  one  source  has  suggested  that 
the  normal  distribution  may  be  considered  as  the  demand  dis- 
tribution at  the  factory  level  and  the  exponential  distribu- 
tion may  be  considered  as  the  demand  distribution  at  the 
wholesale  and  retail  levels,  and  so  on.   Five  characterized 
forms  of  demand  distributions  including  symmetric,  right-skewed. 
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and  left-skewed  will  be  selected  to  test  the  performance  of 
the  candidate  estimators  for  generality,  and  are  shown  in 
Figure  6. 


Beta 

Uniform 
Normal 
Gamma 
Exponential 


0  Demand ,  X 

Figure  6.   Characterized  Forms  of  Demand  Distributions 

Demand  distributions  used  in  the  simulation  and  their 
parameters  are  shown  in  Table  I .   The  parameters  for  the 
normal  distribution  were  selected  so  that  negative  demand 
should  not  be  generated. 

Table  I.   Distributions  and  Their  Parameters 


Distribution 

Parameters 

Uniform 

The  Range  over  [0,30] 

Normal 

Mean  35  and  Variance  100 

Exponential 

Mean  2  0 

Gamma 

a  =  2  and  A  =  0 . 1 

Beta 

A  =  15  and  B  =  5 
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For  each  demand  distribution,  the  99.9th  percentile  of 

demand  except  uniform  distribution  will  be  adopted  as  D 

^        max 

to  apply  the  minimax  rule  under  uncertainty.   Quantile  values 
q  =  ^^/(Cg  +  C^)  will  be  0.1,  0.3,  0.5,  0.7,  and  0.9  for  each 
distribution. 

C.   ALGORITHM  FOR  SIMULATION 

The  program  steps  are  ordered  in  the  following  way: 

1.  Define  a  new  demand  distribution. 

2.  Assign  the  quantile  value  and  maximum  demand. 

3.  Compute  the  optimum  in  cases  of  uncertainty  and  risk. 

4.  Generate  three  random  deviates  from  the  specified 
demand  distribution,  representing  demand  for  the  first 
three  periods. 

5.  Obtain  the  order  statistics. 

6.  Compute  the  estimates  using  the  candidate  estimators. 

7.  Generate  a  random  deviate  for  the  specified  demand 
distribution  for  the  next  period,  compute  the  cost  for 
the  period  and  save  it. 

8.  If  the  period  is  beyond  the  final  period,  go  to  Step  9; 
otherwise  go  to  Step  5. 

9.  If  the  number  of  replications  is  greater  than  the  maxi- 
mum, go  to  Step  10;  Otherwise  go  to  Step  4. 

10.  Find  the  average  cost  for  each  period  over  total 
replications. 

11.  Compute  the  average  cumulative  cost,  the  mixed  cost, 
and  their  relative  cost  deviation  to  the  risk  cost. 
Then  tabulate  the  results. 
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12.  If  quantile  value  is  greater  than  0.9,  go  to  Step  13; 
otherwise  go  to  Step  2. 

13.  If  we  want  to  try  another  demand  distribution  go  to 
Step  1;  otherwise  stop. 

D.   DISCUSSION  OF  THE  SIMULATION  RESULTS 

In  this  section,  the  simulation  results  will  be  displayed 
and  evaluated.   As  mentioned  previously,  the  simulation  run 
was  performed  over  50  replications,  each  of  which  consisted 
of  50  consecutive  demand  periods.   From  these  trials,  the 
average  cumulative  costs  and  the  mixed  costs  were  computed 
for  various  estimators. 

Table  II  shows  the  average  cumulative  cost  for  50  periods 
for  the  given  quantile  values  and  specified  demand  distribu- 
tion.  Errors  in  estimating  D    were  not  considered.   This 

^   max 

performance  of  quantile  estimation  may  be  compared  directly 
with  "best  case"  minimax  results.   Within  a  distribution  we 
may  conclude  one  estimator  is  better  than  the  other  but  it 
is  not  easy  to  conclude  for  all  distributions  using  Table  II. 

Now,  a  measure  of  effectiveness  (MOE)  to  evaluate  the 
performance  of  the  estimators  for  all  distributions  is  needed. 
Let  us  define  the  relative  cost  deviation  of  an  estimator  as 
follows: 


Relative  cost  deviation   _   (cost  using  estimator)  -  (risk  cost) 
(to  the  risk  cost)  (risk  cost) 
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Table  II 

Average  Cumulative  Costs  for  50 
Periods  and  50  Replications 


N^esc 

DisfiK 

Estimator 

0.1 

0.3 

Quantile 
0.5 

(q) 

0.7 

0.9 

Risk 

5.498 

12.737 

15.038 

12.456 

5.310 

Uniform 
(0,3) 

Minimax 
X 

.q 

X 

^q 

X 

_q 

5.498 
5.778 

12.737 
13.624 

15.038 
16.468 

12.456 
13.833 

5.310 
6.290 

5.846 

14.095 

17.397 

15.117 

7.702 

5.782 

13.369 

15.935 

13.174 

7.253 

5.672 

13.366 

16.195 

13.541 

7.011 

X 

q 

5.687 

13.376 

16.205 

13.580 

6.604 

Risk 

7.128 

14.015 

15.997 

13. 805 

6.937 

Normal 
(35,100) 

Minimax 
X 

.q 

X 

::q 
X 

_q 

X 

_q 

X 

q 

11.511 
8.272 

19.943 
15.686 

16.648 
17.818 

15.399 
15.120 

9.426 
8.006 

8.790 

16.870 

19.473 

16.892 

9.417 

7.755 

14.724 

16.868 

14.577 

9.474 

7.934 

15.288 

17.586 

15.010 

8.835 

7.944 

15.232 

17.502 

14.980 

8.329 

Risk 

7.444 

19.657 

27.541 

28.892 

18.642 

Exponen- 
tial 

(Mean  = 
20) 

Minimax 
X 

.q 

X 
.q 

X 

_q 

!^q 
X 

q 

11.705 
7.958 
7.949 

42.088 
20.447 
20.425 

60.824 
28.589 
28.853 

57.014 
30.453 
30.434 

26.512 
20.906 
21.114 

8.389 

21.623 

29.833 

33.840 

25.094 

8.045 

20.576 

28.591 

30.746 

20.950 

8.046 

20.568 

28.547 

30.447 

20.072 

Risk 

6.714 

16.287 

21.141 

20.748 

12.496 

Gamma 
(a  =  2 
A  =  0.1) 

Minimax 

^q 
X 

=q 

X 

_q 

2q 

7.024 
7.253 
7.303 

23.094 
16.963 
17.180 

36.970 
22.275 
22.496 

37.240 
22.141 
22.442 

18.098 
13.649 
14.138 

7.641 

17.761 

22.861 

23.690 

15.932 

7.305 
7.297 

17.067 
17.043 

22.100 
22.098 

21.992 
21.964 

13.863 
13.299 
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Table   II    (Continued) 


>v^pesc 

Quantile 

(q) 

Estimator 

0.1 

0.3 

0.5 

0.7 

0.9 

Risk 

0.071 

0.140 

0.154 

0.129 

0.059 

Beta 

Minimax 

0.260 

0.545 

0.520 

0.246 

0.059 

(A=  15 
B=  5) 

X 

0.079 

0.147 

0.162 

0.135 

0.064 

X 

:q 

X 

_q 

X 

=q 

0.079 

0.150 

0.167 

0.141 

0.070 

0.082 

0.150 

0.163 

0.136 

0.082 

0.078 

0.147 

0.161 

0.135 

0.071 

0.078 

0.147 

0.161 

0.135 

0.067 
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This  relative  cost  deviation  may  be  used  as  an  MOE  to  evalu- 
ate the  performance  of  an  estimator.   The  best  performance 
was  given  by  the  risk  solution  during  the  simulation,  an 
expected  result  with  regard  to  theoretical  considerations. 

Table  III  exhibits  the  relative  cost  deviation  for 
quantile  values  from  0.1  to  0.9  for  each  demand  distribution 
investigated.   As  a  special  case,  if  we  pursue  the  uniform 
demand  distribution  in  Tables  II  and  III,  the  risk  and 
minimax  approach  would  match  each  other  because  the  risk 

cost  solution  X   such  that  F(X  )  =  C_/(C  + C  )  is  the  same 

q  q     o^   s    o 

as  the  minimax  cost  solution  given  by  [C  /(C  +C  )]D 

^      -^    o'   s   o'  •■  max 

whenever  the  maximum  demand  D    is  determined  exactly.   If 

max  "^ 

we  read  through  the  third  column  of  overall  distribution  in 
Table  III,  an  estimator,  X   shows  the  best  performance  among 
candidate  estimators  at  quantile  0.1.   Also,  X  did  better 
at  quantile  0.9.   Excluding  these  two  points,  the  estimator 
X   gave  the  best  performance  at  all  other  quantiles.   Within 
the  uniform  distribution  the  minimax  approach  shows  the  best 
performance  (since  it  is  identical  to  the  risk  solution), 
but  it  shows  the  worst  performance  in  the  other  distribution. 
Similarly,  we  can  say  which  estimator  represents  the  best 
performance  within  a  distribution  by  reading  through  row 
by  row,  and  which  estimator  shows  the  best  performance  at 
a  quantile  by  reading  through  column  by  column. 

We  are  interested  in  the  estimator  which  gives  us  the 
best  performance  for  all  quantiles  over  all  distributions. 
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Table  III 
Relative  Cost  Deviation  From  The  Average  Cumulative  Cost 


"N^Desc. 

Dist'nX 

Estimator 

Quantile  (q) 

Overall 

0.1 

0.3 

0.5 

0.7 

0.9 

Quantile 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Uniform 

Minimax 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

(0.30) 

X 
X 

_q 

X 

_q 

0.0509 

0.0696 

0.0949 

0.1105 

0.1845 

0.5104 

0.0634 

0.1066 

0.1569 

0.2136 

0.4505 

0.9909 

0.0517 

0.0496 

0.0597 

0.0576 

0.3658 

0.5844 

0.0318 

0.0493 

0.0770 

0.0871 

0.3203 

0.5655 

0.0345 

0.0502 

0.0776 

0.0902 

0.2435 

0.4960 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Normal 

MiniiTQX 

0.6149 

0.4229 

0.0405 

0.1154 

0.3589 

1.5526 

(35,100) 

X 

.q 

X 

~q 

X 

_q 

X 

_q 

X 

0.1605 

0.1192 

0.1136 

0.1962 

0.1542 

0.6427 

0.2331 

0.2037 

0.2173 

0.2236 

0.3576 

1.2353 

0.0879 

0.0505 

0.0544 

0.0559 

0.3658 

0.6145 

0,1130 

0.0908 

0.0994 

0.0782 

0.2736 

0.6640 

0.1145 

0.0868 

0.0941 

0.0850 

0.2008 

0.5812 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

EjqxDnen- 

Minimax 

0.5724 

1.1411 

1.2085 

0.9733 

0.4222 

4.3175 

tial 
(Mean  = 
20) 

X 

=q 

X 

_q 

X 

=q 

X 

q 

0.0690 
0.0678 

0.0402 
0.0391 

0.0381 
0.0476 

0.0540 
0.0603 

0.1215 
0.1326 

0.3228 
0.3474 

0.1269 

0.1000 

0.0832 

0.1713 

0.3461 

0.8275 

0.0806 

0.0467 

0.0381 

0.0642 

0.1238 

0.3534 

0.0809 

0.0463 

0.0365 

0.0538 

0.0767 

0.2942 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Gamma 

Minimax 

0.0412 

0.4180 

0.7488 

0.7949 

0.4483 

2.4562 

(a  =  2 
X   =  0.1) 

X 

q 

X 

~q 

X 

q 

X 

_q 

X 

q 

0.0803 
0.0878 

0.0415 
0.0548 

0.0536 
0.0641 

0.0671 
0.0817 

0.0922 
0.1314 

0.3347 
0.4198 

0.1381 

0.0905 

0.0814 

0.1418 

0.2749 

0.7267 

0.0881 

0.0479 

0.0454 

0.10600 

0.1094 

0.3508 

' 

0.0869 

0.0464 

0.0453 

0.0586 

0.0642 

0.3014 
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Table   III    (Continued) 


'^vDesc. 

Dist'n""\ 

Estimator 

Quantile  (q) 

Overall 

0.1 

0.3 

0.5 

0.7 

0.9 

Quantile 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Beta 

Minimax 

2.6322 

2.9002 

2.3714 

0.9108 

0.0085 

8.8291 

(A  =15 
B=  5) 

X 

0.1023 

0.0519 

0.0549 

0.0481 

0.0864 

0.3436 

X 

X 

_q 

X 

q 

0.1123 

0.0753 

0.0877 

0.0952 

0.1823 

0.5528 

0.1481 

0.0736 

0.0558 

0.0575 

0.3852 

0.7202 

0.0916 

0.0508 

0.0459 

0.0468 

0.2019 

0.4368 

\ 

0.0917 

0.0481 

0.0462 

0.0451 

0.1307 

0.3618 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Overall 

Minimax 

3.8657 

4.8822 

4.3752 

2.7944 

1.2379 

17.1554 

Dist'ns 

^q 

0.4630 

3.3224 

0.3551 

0.3749 

0.6388 

2.1542 

X 

~q 

X 

_q 

X 

_q 

X 

q 

0.5644 

0.4795 

0.5736 

0.6744 

1.2543 

3.5462 

0.5527 

0.3642 

0.3345 

0.4841 

1.7378 

3.4733 

0.4049 

0.2855 

0.3058 

0.3453 

1.0290 

2.3705 

0.4085 

0.2771 

0.2997 

0.3327 

0.7159 

2.0346 
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One  way  to  compare  is  to  note  that  the  performance  of  an 

estimator   overall  quantiles  over  all  distributions  can  be 

represented  as  the  total  siim  of  the  relative  cost  deviations. 

Looking  at  the  last  column  for  overall  quantiles  over  all 

distributions  it  may  be  concluded  that  the  best  performance 

(minimum  sum)  is  given  by  the  estimator  X  =  ^[X,   ,.  +X,  > 

q    3   (r-1)    (r) 

"^  ^(    +l^^  ^°^  ^-'-■^  quantiles  over  all  distributions  as  shown 
in  Table  III,  because  the  minimum  relative  cost  deviation  is 
the  closest  to  the  risk  cost  (ideal  cost  solution) .   For 
specific  quantiles,  of  course,  occasionally  other  estimators 
performed  better.   In  an  application  of  the  newsboy  problem 
under  initial  uncertainty,  one  has  information  about  the 
quantile  q  but  often  not  about  the  form  of  the  distribution. 
Accordingly,  the  information  at  the  foot  of  Table  III  may 
provide  information  in  selecting  an  estimator. 

E.   DISCUSSION  OF  SWITCHING  RULES 

When  solving  a  newsboy  problem  repeatedly  under  uncer- 
tainty conditions  we  should  probably  use  the  minimax  rule  for 
the  first  period  since  there  is  no  applicable  amount  of  data 
to  use  in  quantile  estimation.   Each  successive  decision 
period  provides  an  additional  piece  of  demand  data  which 
will  contribute  to  quantile  estimation.   Therefore  we  are 
interested  in  the  number  of  decision  periods  that  should 
elapse  before  we  apply  the  proposed  estimator  X   in  place  of 

Si 

the  minimax  rule.   From  the  simulation  results.  Table  IV 
shows  the  mixed  cost  with  switching  at  period  25,  the  middle 
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Table  IV 
The  Mixed  Cost  for  the  25-th  Period 


Dist>\ 

Estimator 

Qiiantile   (q) 
0.1             0.3             0.5             0.7             0.9 

Unifonn 
(0.30) 

Risk 

Minimax 

X 
X 

_q 

X 

f 

q 

274.879     636.863     751.893     622.810     265.514 
274.879     636.863     751.893     622.810     265.514 
280.729     643.679     761.937     634.742     272.186 
282.951     648.889     766.665     640.032     277.903 
280.185     641.266     759.937     631.729     269.095 
279.242     643.131     760.970     633.978     271.334 
279.595     643.109     761.187     634.082     271.463 

Normal 
(35,100) 

Risk 

Miniinax 

X 

q 

X 

~q 

X 

_q 

X 

-q 

X 

356.404     700.766     799.844     690.275     346.833 
575.550     997.144     832.229     769.934     471.301 
461.583     828.894     820.863     738.497     413.906 
477.707     835.620     825.044     742.866     420.868 
455.163     825.489     819.134     736.412     410.975 
460.505     828.510     820.079     737.839     412.314 
460.535     828.517     820.258     737.909     412.595 

Exponen- 
tial 

(Mean  = 
20) 

Risk 

Minimax 

X 

.q 

X 

=q 

X 

_q 

372.210     982.875  1377.038  1444.599     932.085 
585.271  2104.407  3041.196  2850.687  1325.620 
458.956  1453.137  2073.317  2028.197  1118.620 
458.697  1453.506  2075.067  2027.887  1115.982 
461.103  1456.723  2076.325  2036.363  1146.552 
458.524  1453.449  2072.446  2027.401  1121.038 
458.556  1453.220  2072.413  2027.286  1119.067 

Gairma 
(a=2 
A  =  0.1) 

Risk 

Miniinax 
X 

^q 

X 

~q 

X 

_q 

X 

q 

335.684     814.349  1057.038  1037.387     624.824 
351.179  1154.714  1848.499  1861.989     904.908 
344.370     983.554  1414.183  1400.872     759.497 
344.799     985.234  1415.761  1402.143     759.026 
348.497     983.682  1410.905  1406.584     768.769 
344.654     982.802  1410.740  1400.742     754.786 
344.504     982.942  1411.690  1400.479     755.480 
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Table 

IV    (Continued) 

\Desc. 

Dist^\ 

Estimator 

Quantile 

(q) 

0.1 

0.3 

0.5 

0.7 

0.9 

Risk 

3.573 

6.993 

7.699 

6.437 

2.948 

Beta 

Minimax 

12.988 

27.275 

26.003 

12.299 

2.973 

(A=15 
5=    5) 

X 

q 

7.691 

15.796 

15.779 

9.190 

3.021 

X 

=q 

X 

_q 

X 

q 

7.726 

15.789 

15.779 

9.250 

3.056 

7.723 

15.840 

15.809 

9.214 

3.015 

7.680 

15.788 

15.781 

9.191 

3.014 

7.678 

15.789 

15.779 

9.190 

3.015 
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period  in  our   fifty-period  simulation.   The  mixed  costs 
here  are  the  sum  of  costs  for  the  first  24  periods  given  by 
the  minimax  rule  and  the  costs  given  by  an  estimator  from 
the  25-th  to  the  last  period.   When  we  examine  various  dis- 
tributions in  Table  IV  the  minimum  mixed  cost  is  obtained 
by  the  risk  approach,  which  is  to  be  expected,  and  the 
proposed  estimator,  X  ,  appears  to  provide  the  minimum 

Si 

mixed  cost  among  candidate  estimators,  although  by  a  narrow 
margin. 

Now,  we  wish  to  find  the  number  of  periods  to  retain 
the  minimax  rule  before  applying  the  proposed  estimator  X  . 
Mixed  costs  were  computed  for  various  switching  periods  to 
find  the  period  which  gives  the  minimum  mixed  cost  for  each 
quantile  of  every  distribution  except  uniform  distribution, 
since  the  proposed  estimator  X  doesn't  give  better  results 
than  the  minimax  rule  for  the  uniform  distribution.   Table 
V  shows  the  best  switching  periods  (minimax  to  estimator  X  ) 
for  the  quantiles  and  distributions  used  in  the  simulation 
during  4  0  periods. 

Table  V 
Best  Switching  Periods  Using  the  Proposed  Estimator  X 


Distribution 

Quantile 

(q) 

0.1 

0.3 

0.5 

0.7 

0.9 

Normal 

4 

7 

27 

11 

11 

Exponential 

6 

4 

4 

4 

4 

Gamma 

13 

5 

4 

4 

6 

Beta 

4 

4 

4 

4 

48 
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Looking  at  Table  V,  we  see  that  the  best  switching  period 
is  27  at  quantile  0.5  of  the  normal  distribution.   This 
result  is  caused  by  the  fact  that  the  minimax  rule  behaves 
well  around  the  median  when  demand  has  a  symmetric  distri- 
bution.  The  occurrence  of  switching  period  4  8  at  quantile 
0.9  of  the  beta  distribution  can  be  explained  by  the  fact 
that  with  the  beta  distribution  used  in  the  simulation, 
the  optimal  S*  obtained  by  the  minimax  rule  is  close  to 
the  ideal  S*  at  this  quantile,  i.e.,  the  minimax  S*  at  0.9 
is  equal  to  0.8767  and  the  ideal  S*  is  0.8660. 

We  shall  suggest  6  periods,  4,  5,  6,  7,  11,  and  13  as 
the  candidates  for  the  switching  period  for  all  quantiles 
over  all  distributions.   Later  periods  are  not  included  in 
the  candidates  since  the  estimator  X   shows  much  better 

q 

performance  (with  reference  to  the  foot  of  Table  III  for 
nearly  all  quantiles  over  all  distributions) .   Table  IV 
shows  the  mixed  costs  obtained  by  the  estimator  X   for  the 

q 

25th  period.   Similar  tables  (not  shown)  were  computed  for 
various  other  periods,  and  showed  that  the  mixed  costs 
increased  as  the  period  increased.   We  propose  the  relative 
mixed  cost  deviation  for  a  candidate  period  using  the  esti- 
mator X   as  an  MOE .   It  may  be  reasonable  to  conclude  that 

q 

the  period  giving  the  minimum  is  the  most  preferable.   Table 
VI  shows  the  relative  mixed  cost  deviation  at  each  candidate 
period  for  all  quantiles  over  all  distributions  using  the 
estimator  X  . 

q 
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Table  VI 
Relative  Mixed  Cost  Deviation  Using  the  Estimator  X 


Switching 
Period 

Distribution 

Estimator 

0.1 

0.3 

Quantile  (q 
0.5    0.7 

) 
0.9    Total 

Uniform 

.0345 

.0502 

.0776 

.0902 

.2435  .4960 

Normal 

.1145 

.0868 

.0941 

.0880 

.2008  .5812 

4 

Exponential 

X 

q 

.0809 

.0463 

.0365 

.0538 

.0767  .2942 

Gamma 

.0869 

.0464 

.0453 

.0586 

.0642  .3014 

Beta 

.0917 

.0481 

.0462 

.0451 

1307  .3618 

Overall 

.4085 

.2778 

.3327 

.3327 

.7159  2.0346 

Uiiform 

.0324 

.0471 

.0675 

.0762 

.1720  .3952 

Nbimal 

.1201 

.0868 

.0751 

.0673 

.1462  .4955 

5 

Exponential 

X 

q 

,0708 

,0548 

.0561 

.0669 

.0851  .3337 

Gamma 

.0576 

.0422 

.0555 

.0713 

.0708  .2971 

Beta 

.1235 

.1003 

.0902 

.0570 

.0868  .4578 

Overall 

.4044 

.3312 

.3444 

.3387 

.5606  1.9793 

Uniform 

.0317 

.0457 

.0626 

.0693 

.1229  .3322 

Normal 

.1255 

.0863 

.0663 

.0644 

.1145  .4570 

6 

Exponential 

1 
q 

.0657 

.0697 

.0742 

.0765 

.0894   .3755 

Gamma 

.0491 

.0427 

.0675 

.0814 

.0774  .3181 

Beta 

.1715 

.1591 

.1400 

.0790 

.0608  .6104 

Overall 

.4435 

.4035 

.4106 

.3706 

.4650  2.0932 

Lftiiform 

.0300 

.0420 

.0547 

.0601 

.1020  .2888 

Normal 

.1311 

.0846 

.0533 

.0572 

.1084  .4346 

7 

Exponential 

X 

q 

.0738 

.0924 

.0980 

.0944 

.0987  .4573 

Gaitita 

.0423 

.0487 

.0754 

.0924 

.0847  .3435 

Beta 

.2306 

.2237 

.1935 

.0973 

.0481  .7932 

Overall 

.5078 

.4914 

.4749 

.4014 

.4419  2.3174 
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Table  VI    (Continued) 


Switching 
Period 

Distribution 

EstiiiHtor 

0.1 

0.3 

Quantlle  (q) 
0.5    0.7 

0.9 

Itotal 

Uniform 

.0262 

.0323 

.0384 

.0422 

.0502 

.1891 

Nbnrel 

.1727 

.0952 

.0368 

.0526 

.0988 

.4561 

11 

Exponential 

X 

q 

.0859 

.1641 

.1779 

.1594 

.1292 

.7165 

Gamma 

.0194 

.0797 

.1253 

.1384 

.1041 

.4669 

Beta 

.4119 

.4456 

.3793 

.1705 

.0375 

1.4448 

Overall 

.7161 

.8169 

.7577 

.5631 

.4196 

3.2734 

Uniform 

.0247 

.0280 

.0323 

.0383 

.0446 

.1679 

Normal 

.1916 

.1047 

.0322 

.0544 

.1070 

.4899 

13 

Exponential 

^q 

.1057 

.2066 

.2228 

.1932 

.1413 

.8696 

Gainra 

.0910 

.0930 

.1502 

.1627 

.1158 

.6127 

Beta 

.5236 

.5666 

.4823 

.2166 

.0362 

1.8253 

Overall 

.9366 

.9989 

.9198 

.6652 

.4449 

3.9654 
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The  relative  mixed  cost  deviation  at  each  candidate 
period  is  shown  in  Table  VII.   The  minimum  deviation  is 
given  at  period  five  and  the  deviation  is  increasing  as  the 
period  increases  from  the  fifth  period.   Using  this  result 
we  may  propose  that  the  switching  period  from  the  minimax 
rule  to  applying  the  estimator  X   should  be  the  fifth  period, 


Table  VII 

The  Relative  Mixed  Cost  Deviation 
At  Each  Candidate  Period 


Period 

Cost  Deviation 

4 

2.0346 

5 

1.9  79  3   (minimum) 

6 

2.0932 

7 

2.3174 

11 

3.2734 

13 

3.9654 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

The  minimax  approach  is  a  well  known  approach  to  decision 
making  under  uncertainty.   For  a  newsboy-type  problem,  we 
developed  five  candidate  estimators  which  could  be  used  as 
data  aggregates,  and  evaluated  them  by  comparing  with  the 
minimax  rule.   These  candidate  estimators  were  tested  for 
quantiles,  0.1,  0.3,  0.5,  0.7,  and  0.9  of  several  distribu- 
tions:  uniform,  normal  (symmetric),  exponential,  gamma 
(right-skewed)  and  beta  (left-skewed) .   Since  all  possible 
quantiles  and  other  kinds  of  distributions  were  not  covered, 
there  may  be  some  amount  of  loss  of  generality. 

With  these  limitations,  as  discussed  in  Chapter  IV,  the 

estimator  X  =  rr-ix ,      ,,+X,  ,+X,  ,,,]  developed  on  the  basis 
q    3   (r-1)    (r)    (r+1)         ^ 

of  order  statistics  showed  the  best  performance  among  those 
considered.   The  simulation  results  suggest  that  if  a  single 
decision  procedure  over  all  quantiles  is  needed,  the  esti- 
mator X   is  recommended  and  the  time  for  switching  to  this 

q 

estimator  should  be  at  the  fifth  period.   The  proposed 
procedure  may  be  summarized  in  the  following  way: 

(1)  Estimate  D 

max 

(2)  Compute  the  quantile  q  =  C^/(Cg  +  C^). 

(3)  Apply  the  minimax  rule  S*  =  q-D_^^  for  the  first 

cc   J.  J.     max 

four  periods. 

(4)  Compute  the  value  r  such  that  r  =  [nq  +  0.5]  from  the 
fifth  period  where  n  is  the  period  at  which  we  wish 
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to  estimate  the  demand  and  [X]  denotes  the  largest 

integer  of  X. 

(5)   Compute  the  order  statistics  X,   ,>,  X,  .,  and  X,  ,,> 
^  (r-1)    (r)        (r+1) 

and  apply  the  proposed  estimator  X  : 

^q  =   3^^(r-l)  "^^(r)  "^^(r+l)^  * 

In  this  study,  it  was  implicitly  assumed  that  the  value 
of  demand  for  each  period  was  obtainable  even  if  D  >  S.   For 
further  study,  the  case  when  demand  data  greater  than  the 
quantity  S  is  not  available  is  worthy  of  investigation. 
Also,  a  non-parametric  approach  which  can  handle  extreme 
quantiles  (e.g.,  0.01  or  0.99)  and  extremely  skewed  distri- 
butions could  be  investigated.   The  comparison  between  a 
parametric  approach  when  the  demand  data  are  assumed  to  fit 
a  probability  distribution  and  non-parametric  approach  is 
also  recommended.   Finally,  future  simulation  work  on  this 
problem  should  also  record  cost  variances,  so  that  a  statis- 
tical  analysis  could  be  made  of  the  simulation  results. 
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APPENDIX 
PROGRAM   LISTING    OF    THE    SIMULATION 
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